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Abstract 

By using the finite temperature quantum field theory, we calculate the finite temperature effective 
potential and extend the improved quark mass density-dependent model to finite temperature. It 
is shown that this model can not only describe the saturation properties of nuclear matter, but also 
explain the quark deconfinement phase transition successfully. The critical temperature is given 
and the effect of uj- meson is addressed. 
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I. INTRODUCTION 



Quark- meson coupling(QMC) model suggested by Guichon [1] is a famous hybrid quark 
meson model, which can describe the saturation properties of nuclear matter and many 
other properties of nuclei successfully. In this model, the nuclear system was suggested as 
a collection of MIT bag, uj- meson and a- meson and the interactions between quarks and 
mesons arc limited within the MIT bag regions because the quark cannot escape from the 
MIT bag. This model had been extended by many authors, for example, including p meson 
to discuss the neutron matter [2] , adding hyperons to study the strange hadronic matter [3] , 
suggesting the density- dependent vacuum energy to investigate the effect of environment 
[4] and etc. [5]. 

Although the QMC model is successful in describing many physical properties of nuclear 
systems, but as was pointed in our previous paper [6], it has two major shortcomings: (1) 
It is a permanent quark confinement model because the MIT bag boundary condition can- 
not be destroyed by temperature and density. It cannot describe the quark deconfinement 
phase transition. (2) It is difficult to do nuclear many-body calculation beyond mean field 
approximation(MFA) by means of QMC model, because we cannot find the free propaga- 
tors of quarks and mesons easily. The reason is that the interactions between quarks and 
mesons are limited within the bag regions, the multirefiection of quarks and mesons by 
MIT bag boundary must be taken into account for getting the free propagators. These two 
shortcomings actually inherited from the MIT bag. 

To overcome these two shortcomings, in our previous paper [6], we suggested an improved 

quark mass density- dependent (QMDD) model. The QMDD model was suggested by Fowler, 

Raha and Weiner [7] many years ago. According to the QMDD model, the masses of u, d, 

s quarks and corresponding antiquarks are given by 

B - 
ruq^- — {i = u,d,u,d) (1) 
3ns 

B 

ms,s = + - — (2) 

where nisQ is the current mass of the strange quark, B is the bag constant, ub is the baryon 
number density 

riB ^\^{nu + nd + ns), (3) 



Uu, rid, IT'S represent the density of u quark, d quark, and s quark, respectively. As was ex- 
plained and proved in Ref. [8, 9], the basic hypothesis (1), (2) correspond to a confinement 
mechanism of quarks because when volume V ^ oo, ub ^ 0, and then rUg oo, quark will 
be confined. Using QMDD model, many authors investigated the dynamical and thermo- 
dynamical properties of strange quark matter and strange star and found that the results 
given by QMDD model are nearly the same as those obtained in MIT bag model [8-16]. But 
a nice advantage for QMDD model is that the MIT bag boundary condition has been given 
up. The ansatz Eq. (1) is employed to replace the MIT bag boundary. 

But QMDD model still has two shortcomings: (1) It is still an ideal quark gas model. No 
interactions between quarks exist except a confinement ansatz Eqs. (1-2). (2) It still cannot 
explain the quark deconfincment phase transition and give us a correct phase diagram as 
that given by lattice QCD because when ub 0,mq — > oo, T — > oo [8]. To overcome the 
difficulty (2), we introduced a new ansatz that nig is not only a function of density, but also 
depends on temperature T [8,12,15], we suggested 



and extended the QMDD model to a quark mass density- and temperature- dependent 
model(QMDTD) model. The ansatz(4) guarantees that — > when T — > Tc- It changes 
the permanent confinement mechanism(MIT bag) to a nonpermanent confinement mecha- 
nism (Friedberg-Lcc(FL) soliton bag) in the QMDTD model. Since the vacuum density B 
equals to the different value between the local false vacuum minimum and the absolute real 
vacuum minimum of nonlinear scalar field in FL model, we introduced a nonlinear scalar 
field to improve the QMDD model in Refs. [17, 22] and changed the ad-hoc ansatz(4) to an 
output B(T) curve which can be calculated from the improved QMDD model. 

To improve the shortcoming(l) of QMDD model, we introduced the cu— meson and a— 
meson in QMDD model to mimic the repulsive and attractive interactions between quarks 
in Refs. [6, 23]. The interaction between quarks and nonlinear a field forms a FL soliton bag 
[17, 22]. The qqcu and qqa interaction guarantees that we can get valid saturation properties, 
the equation of state and the compressibility of nuclear matter [6]. But the deconfincment 
phase transition has not yet been studied in Ref. [6]. This motivates us to study the 
deconfincment properties of the improved QMDD(IQMDD) model in this paper. We would 



B = Bo[l - {T/Tcfl 0<T<Tc, 



(4) 
(5) 



:T>Tc 
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like to emphasize that this is the basic important advantage for IQMDD model, because 
the saturation properties can be explained not only by IQMDD model but also by QMC 
model. The reason for the explaination of quark deconfinement by IQMDD model is that 
MIT boundary constraint has been dropped and interactions between quark and mesons 
have been extended to the whole space. Since instead of the MIT bag in QMC model, a 
FL soliton bag is introduced in IQMDD model, we can use our model to discuss the quark 
deconfinement phase transition. The spontaneous breaking symmetry of nonlinear a field 
will be restored and the soliton bag will disappear at critical temperature [24-27]. 

The organization of this paper is as follows. In the next section, we give the main formulae 
of the IQMDD model and effective potential at finite temperature. The soliton solutions 
of IQMDD model at different temperature and other numerical results are presented in the 
third section. The last section contains a summary and discussions. 



II. THE IQMDD MODEL AT ZERO AND FINITE TEMPERATURE 



The effective Lagrangian density of the IQMDD model is given by 



C = ipli-f^'df, -m^- fa- gj''ujf,]ip 



(6) 



where F^i, = d^uj^ — d^uj^, tp represents the quark field, iriq = is the mass of u{d) quark, 
the a and u field are not dependent on time, / is the coupling constant between the quark 
field ip and the scalar meson field a, g is the coupling constant between the quark field ip 
and the vector meson field cj^, U{a) is the self interaction potential for a field. We omit 
the contribution of the s quark and consider the nuclear system only in this paper. The 
potential field f/(cr) is chosen as [18] 



^ ' 2! 3! 

6^ > 3ac 



4! 



(7) 
(8) 



The condition ([8]) ensures that the absolute minimum of U{a) is at a = a,^ ^ 0. The 
potential U (a) has two minima: one is the absolute minimum 
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it corresponds to the physical vacuum, and the other is at ctq = 0, it represents a metastable 
local false vacuum. We take U{av) — and the bag constant B can be expressed as 

From Eq. (6), we obtain the equation of motion for quark as 

{lYd^ -m,-fa- 57^0;^)^ = 0, (11) 
and the equations for the scalar meson field and vector meson field as 

d,d^a+^^-m, (12) 

d.F"^ + m^LU^ = giPj^ip. (13) 

respectively. Using an approximation as that of the QMC model, we replace a{r, t) — > cr(r), 
u!fi{r,t) SfiOUj{r) and consider a fixed occupation number of valence quarks (3 quarks for 
nucleons, and quark- ant iquark pair for mesons)only. In the following, we will discuss the 
ground state solution of the system. The Hamiltonian density is 

7i = • V + Pirug + fa) + + \K 

+ 1 (Va)^ + - \{^u^f - \mlu,u'^. (14) 

where a. and (3 are the Dirac matrices Ii„ is conjugate field of the scalar meson field. One can 
construct a Fock space of quark states and expand the operator in terms of annihilation and 
creation operators on the space with spinor function c/?^, which satisfies the Dirac equation 
[17]: 

[d-p + (3{mq + fa)+ guj]'^^ = ±en^^. (15) 

The functions satisfies the normalized condition j if^ip^d^r — 1. From Eq. (14), the 
total energy of the system is given by 

E{a,u;)^J2'r^+ /[^(Va)^ + t/(a)-i(Vu;)^-^m>Vr. (16) 

n 

Substituting Eq. (15) into Eqs. (12, 13), and using the variational principle under the 
spherical symmetric condition, we find when a and cu satisfy the follow equations 



da 

n 



-W^u + mlu ^ g^ip^l^ipn- (17) 

n 

respectively where if^ — (fin Jo, we have the minimum of E{a,u!). 

We discuss the ground state solution of the system now. The quark spinor in the lowest 
state is assumed [17-19]: 

= f ^^^^\ \ Xm, Xm= i ]or\^], (18) 

where a are the Pauli matrices. Substituting Eq. (18) into Eq. (15), we get the equations 
of spinor components u and v as 

du{r) 



dr 



= -[e + mg + fa{r) - guj]v{r), 



dv(r) 2 

--v{r) + [e - 777,5 - /^(^) - g^^]u{r). (19) 

respectively. The normalized condition then reads as Att J^[u'^{r) + v'^{r)]r'^dr — 1. 

From Eq. (17) and Eq. (18), after summation of the quark states, we obtain the equation 
of motion of the a, ut field 

d'^a 2 da dU(a) 9 os ,^^s 

d^u 2 du 



dr^+r^- <^ = -^9{u' + ^ ) ^ ^(^)- (21) 
where, the number of quarks is A?" = 3 for baryons and A?" = 2 for mesons. In the following 
discussions, we only consider the case N — 3. To get a self-consistent solution of Eqs. (19, 
20, 21), we select the boundary conditions for quark field and a field as [17, 22] 

v{r = 0) = 0, u{r — 00) = 0, 
a'{r = 0) = 0, a{r = 00) = a^. (22) 

respectively. Noting that the r ^ 00 asymptotic behavior of the cu field given by Eq. (21) 
tends to an exponent decay wave because F{r 00) — * for a soliton bag, we can find the 
corresponding Green function G{r,r') easily and obtain the lu field by integral [20]: 

PCX) 

uj{r)= r'^dr'F{r')G{r,r') (23) 
Jo 

The numerical results of u{r), v{r), a{r) and a;(r) will be shown in next section. 



In order to study the deconfinement phase transition, we turn to extend IQMDD model 
to finite temperature. The appropriate framework is the finite temperature quantum field 
theory. The finite temperature effective potential plays a central role within this framework. 
Under mean field approximation, the meson field operators can be replaced by their expec- 
tation values, —>■ u!fj_ = Shqlu = 5^o;^Pbo[20, 21]. Using the method of Dolan and Jackiw 
[24], up to one- loop approximation, the effective potential reads: 

V{a- T; fi- V^) = U{a) + T) + V>((7; T; /x; V^), (24) 

where 

V^ = gu}^ -^pBo (25) 

is the contribution of a;- field, pso is saturation density of nuclear matter, T is temperature 
and 



Vb{<j; T)^^j^" dxx'^ln (l - e-V(-^+-l/^^)) , (26) 

T-'4 pco 

Vpia; T; p; V^) = -12 ^ ^ j dxxHn (l + e-(V(-^+-^"/^^)-/^"/^+^-/^)) , (27) 

where the minus sign of Eq. (27) is the consequence of Fermi-Dirac statistics. The degenerate 
factor 12 comes from: 2 (particle and antiparticle), 2 (spin), 3 (color), m^ and are the 
effective masses of the scalar field a and the quark field respectively: 

ml = a + ha{T) + ^^'(T) (28) 

m,^^ + fa{T). (29) 

Ug 

We see from Eqs. (24-27) that the scalarlike interaction fip^^aip gives contribution to 
effective masses of quark and a meson and then forms a confined soliton bag, the vectorlike 
interaction gip'^^^cun'ip gives contribution to an effective chemical potential of quarks. In 
fact, this finite temperature effective potential for FL model had been calculated by many 
others authors including us [22, 25-27] , but without uj- field. 

In the soliton bag model, the finite temperature vacuum energy density B{T) is defined 

as 

B{T; p; V^) = F(ao; T; p; K,) - V{a,; T; p; V^). (30) 



It is the different from the values at the perturbative false vacuum state and the values 
at the physical real vacuum state of the finite temperature effective potential. At critical 
temperature Tq of quark deconfinement phase transition, B equal to zero: B{Tc) = 0. 

III. RESULTS 

Before numerical calculations, we consider the parameters of IQMDD model at first. To 
guarantee our model can be used not only to explain the saturation properties of nuclear 
matter, but also to describe the quark deconfinement phase transition, we fix our parameters 
as those in Ref. [6]. Hereafter we fix the parameters as: at zero temperature, the bag 
constant S=174 MeV/m~^, the masses of uj- meson and a- meson = 783 MeV and 
m„ = 509 MeV, the coupling constant / = 5.45 and g = 3.37 respectively and the mass of 
nucleon Mjv = 939 MeV, b = —8400 MeV. As was shown in Ref. [6], this set of parameters 
can give us reasonable properties of nuclear matter, including the value of saturation point, 
the equation of state and the compressibility, we will prove in this section that we can also 
obtain a reasonable quark deconfinement critical temperature Tq by means of this set of 
parameters. 

The set of coupled differential equations (19), (20) with boundary conditions Eq. (22) 
can be solved numerically at zero temperature. Our results at zero temperature are shown 
in Figs. (1-4). The variation of the scalar field o" as a function of the radius r is presented in 
Fig. 1. We see that the value of a inside the hadron is very different from that of outside: 
inside a is less than zero, and outside cr — > cr^. The transitional values of a field through 
the surface from inside to outside is abrupt. The curve in Fig. 1 is very similar as that 
given by Refs. [17, 22] where the uj- meson is omitted. This is reasonable because the main 
contribution to form a soliton bag comes from the nonlinear a field. The variation of wave 
function of quark field is shown in Fig. 2 and Fig. 3. In Fig. 2, we plot the curves of quark 
wave function u, v vs. r respectively. In Fig. 3, we plot the curve of quark density — 
vs. r. The soliton bag is exhibited in this figure transparently. The curve of cu- field vs. r 
given by Eq. (21) is shown in Fig. 4. This solution is obtained from Green function method 
and Eq. (23). We see from this Figure that the curve of a;- field decays considerably when 
r is large than the soliton bag radius. 

Now we turn to investigate the case of finite temperature. The effective potential at finite 
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temperature can be obtained by numerical calculations from the set of Eqs. (24-27). The 
curves of effective potential U{a;T; ijl;V^) at different temperatures T — MeV, T = 80 
MeV, T = 127 MeV and T = 204 MeV are shown in Fig. 5 respectively. The shape of 
the effective potential confirms that a first order phase transition will take place [22, 25-27]. 
We see from Fig. 5 there are two vacua where one corresponds to the physical vacuum and 
the other corresponds to the false vacuum when T < 127 MeV. But when T increases to 
a critical temperature Tc — 127 MeV, these two vacua degenerate and the values of bag 
constant tends to zero. 

To show this behaviour more clearly, we plot the bag constant B vs. T curve by solid 
curve in Fig. 6, we see that the bag constant B decreases as temperature T increases. When 
T approaches to Tc, B approaches to zero and quark dcconfincd phase transition happens. 

When temperature increases to the regions 127 MeV < T < 204 MeV, as shown in Fig. 
5 the physical vacuum becomes instable and the purturbative vacuum becomes stable. The 
soliton solution tends to disappear [22]. When T approached to 204 MeV, the effective 
potential becomes to have an unique minimum only. Such potential no longer ensures the 
existence of the soliton bag anymore because the spontaneous breaking symmetry has been 
restored. Wc have zero solution a = in the regions T > 204 MeV only. 

To illustrate our soliton solutions more transparently, we plot the soliton solutions by 
fixed the temperature T = MeV, T = 80 MeV and T = 125 MeV in Fig. 7. It is shown in 
Fig. 7 that the curves stretch slowly to infinite with increasing temperature. The radius of 
soliton increases and the skin of the soliton becomes thicker when the temperature increases. 
It is of course very reasonable. 

All characters of the soliton solutions found by IQMDD model are almost the same as 
those given by Rcf. [22] where the uj field has not been take into account. But we hope to 
emphasize although the uj field will not affect the main characters of the soliton solutions, but 
it will change their detailed behavior. To address the effect of u field on the deconfinement 
transition, as an example, we plot the B{T) curve without cu field by dotted line in Fig. 6. 
We find the critical temperature from the condition B{Tc) = changes to Tc = 140 MeV 
which is higher than 127 McV where the uo field existed. If we notice that the interaction 
of qqoj is repulsive, and then it can help to deconfine quarks, the decrease of Tc due to uj- 
field is natural. 
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IV. SUMMARY AND DISCUSSION 



By using the finite temperature quantum field theory, we calculate the finite temperature 
effective potential and extend our previous discussions to finite temperature. It is shown 
that the IQMDD model can not only describe the saturation properties of nuclear matter, 
but also explain the quark deconfinement phase transition successfully. We have shown 
the soliton solution curves for different temperatures and found the critical temperature 
Tc — 127 MeV. The uj- field in IQMDD model is important. The repulsive qquj interaction 
plays the central role to describe the saturation properties, and affect the critical temperature 
remarkably. Comparing to naive QMC model, the advantage of IQMDD model is obvious 
because instead of MIT bag, a Priedberg-Lee sohton bag exists in this model. 

Of course, there are still shortcomings in the IQMDD model. For example, the interac- 
tions between quarks and mesons are still isospin independent. The chiral phase transition 
cannot be discussed by this model because it is lack of chiral symmetry. To Overcome these 
shortcomings, we hope to add the p- meson and tt- meson in IQMDD model in the near 
future. 
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FIG. 1: The a field as a functions of r at zero temperature. 
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FIG. 2: Quark wave functions in arbitrary unit as a functions of r. 
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FIG. 3: The quark density «^(r) — v'^{r) in arbitrary unit as a functions of r. 
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FIG. 4: The w field as a functions of r. 
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FIG. 5: The temperature-dependent effective potential. 
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FIG. 6: The bag constant B(T) as functions of T, the sohd hne corresponds to IQMDD model 
with w-meson and the dotted hne corresponds to the same model without w-meson. 
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FIG. 7: The soliton solutions for different temperature T=0 MeV, 80 MeV, 125 MeV 
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